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Linear elastic behavior 


A linear elastic material model: See Also 
e is valid for small elastic strains (normally less than 5%); About the material library 
e can be isotropic, orthotropic, or fully anisotropic; Elastic behavior 
e can have properties that depend on temperature and/or other field variables; and 
e can be defined with a distribution for solid continuum elements in In Other Guides 
Abaqus/Standard. 
*ELASTIC 


The following topics are discussed: 
e Defining linear elastic material behavior 
e Directional dependence of linear elasticity 
e Stability of a linear elastic material 
e Defining isotropic elasticity 
e Defining orthotropic elasticity by specifying the engineering constants 
e Defining transversely isotropic elasticity 
e Defining orthotropic elasticity in plane stress 
e Defining orthotropic elasticity by specifying the terms in the elastic stiffness matrix 
e Defining fully anisotropic elasticity 
e Defining orthotropic elasticity for warping elements 
e Defining elasticity in terms of tractions and separations for cohesive elements 
e Defining isotropic shear elasticity for equations of state in Abaqus/Explicit 
e Defining the elastic transverse shear modulus 
e Elements 


Creating a linear elastic 
material model 


Products: Abaqus/Standard Abaqus/Explicit Abaqus/CAE 


Defining linear elastic material behavior 


The total stress is defined from the total elastic strain as 
el el 
a0 =D"%e", 


where @ is the total stress (“true,” or Cauchy stress in finite-strain problems), D® is the fourth-order elasticity tensor, and et! 
is the total elastic strain (log strain in finite-strain problems). Do not use the linear elastic material definition when the elastic 
strains may become large; use a hyperelastic model instead. Even in finite-strain problems the elastic strains should still be 
small (less than 5%). 


Defining linear elastic response for viscoelastic materials 


The elastic response of a viscoelastic material (Time domain viscoelasticity) can be specified by defining either the 
instantaneous response or the long-term response of the material. 


Instantaneous response 


To define the instantaneous response, experiments to determine the elastic constants have to be performed within time spans 
much shorter than the characteristic relaxation time of the material. 


Input File Usage: 
*ELASTIC, MODULI=INSTANTANEOUS 
Abaqus/CAE Usage: 


Property module: material editor: Mechanical > Elasticity > Elastic: Moduli time scale (for viscoelasticity) : 
Instantaneous 


Long-term response 


If, on the other hand, the long-term elastic response is used, data from experiments have to be collected after time spans 
much longer than the characteristic relaxation time of the viscoelastic material. Long-term elastic response is the default 
elastic material behavior. 


Input File Usage: 
*ELASTIC, MODULI=LONG TERM 


Abaqus/CAE Usage: 
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Property module: material editor: Mechanical > Elasticity > Elastic: Moduli time scale (for viscoelasticity): Long-term 


ry 


Directional dependence of linear elasticity 


Depending on the number of symmetry planes for the elastic properties, a material can be classified as either isotropic (an 
infinite number of symmetry planes passing through every point) or anisotropic (no symmetry planes). Some materials have a 
restricted number of symmetry planes passing through every point; for example, orthotropic materials have two orthogonal 


symmetry planes for the elastic properties. The number of independent components of the elasticity tensor D* depends on 
such symmetry properties. You define the level of anisotropy and method of defining the elastic properties, as described 
below. If the material is anisotropic, a local orientation (Orientations) must be used to define the direction of anisotropy. 


rt 


Stability of a linear elastic material 


Linear elastic materials must satisfy the conditions of material or Drucker stability (see the discussion on material stability in 


Hyperelastic behavior of rubberlike materials). Stability requires that the tensor D*@ be positive definite, which leads to 
certain restrictions on the values of the elastic constants. The stress-strain relations for several different classes of material 
symmetries are given below. The appropriate restrictions on the elastic constants stemming from the stability criterion are 
also given. 


r 


Defining isotropic elasticity 


The simplest form of linear elasticity is the isotropic case, and the stress-strain relationship is given by 


Cae) 1/E -v/E-v/E 0 0 0 _ 
£99 —v/E 1/E -v/E 0 0 0 022 
£33 _v/E-v/E 1/E 0 0 0 oa 
wf | 0 0 0 1/46 0 0 on 
13 0 0 0 0 1/4 0 013 
23 0 0 0 0 oO 1/4 023 


The elastic properties are completely defined by giving the Young's modulus, E, and the Poisson's ratio, ’. The shear modulus, 
G, can be expressed in terms of FE and v as G = E/2(1+ 1). These parameters can be given as functions of temperature 
and of other predefined fields, if necessary. 


In Abaqus/Standard spatially varying isotropic elastic behavior can be defined for homogeneous solid continuum elements by 
using a distribution (Distribution definition). The distribution must include default values for E and v. If a distribution is used, 
no dependencies on temperature and/or field variables for the elastic constants can be defined. 


Input File Usage: 
*ELASTIC, TYPE=ISOTROPIC 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Elasticity > Elastic: Type: Isotropic 
Stability 


The stability criterion requires that # > 0, G > 0, and —1 < v < 0.5. Values of Poisson's ratio approaching 0.5 result in 
nearly incompressible behavior. With the exception of plane stress cases (including membranes and shells) or beams and 
trusses, such values require the use of “hybrid” elements in Abaqus/Standard and generate high frequency noise and result in 
excessively small stable time increments in Abaqus/Explicit. 


In Abaqus/Standard it is recommended that you use solid continuum hybrid elements for linear elastic materials with Poisson's 
ratio greater than 0.495 (i.e., the ratio of K/p greater than 100) to avoid potential convergence problems. Otherwise, the 


analysis preprocessor will issue an error. You can use the “nonhybrid incompressible” diagnostics control to downgrade this 
error to a warning message. 


Input File Usage: 
Use the following option to downgrade an error message to a warning message: 


*DIAGNOSTICS, NONHYBRID INCOMPRESSIBLE=WARNING 


rt 


Defining orthotropic elasticity by specifying the engineering constants 
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Linear elasticity in an orthotropic material is most easily defined by giving the “engineering constants”: the three moduli £1, 
Ey, E3; Poisson's ratios 112, 113, Y23; and the shear moduli G32, G13, and G23 associated with the material's principal 
directions. These moduli define the elastic compliance according to 


E11 1/E,  —V_,/E_ —v3,;/E3 0 0 0 ou 
E22 —V"y2/E,  1/E, —v3/E; 0 0 0 022 
£33 —113/Ey —V23/E_  1/B3 0 0 0 033 
mw ( 0 0 0 1G», 0 0 on 
13 0 0 0 0 1/Gi, 0 013 
23 0 0 0 0 0 1/43 023 


The quantity 14; has the physical interpretation of the Poisson's ratio that characterizes the transverse strain in the j-direction, 
when the material is stressed in the /-direction. In general, 14; is not equal to ¥3;: they are related by V3 / Bg=V 54 / Ej. The 
engineering constants can also be given as functions of temperature and other predefined fields, if necessary. 


In Abaqus/Standard spatially varying orthotropic elastic behavior can be defined for homogeneous solid continuum elements 
by using a distribution (Distribution definition). The distribution must include default values for the elastic moduli and 
Poisson's ratios. If a distribution is used, no dependencies on temperature and/or field variables for the elastic constants can 
be defined. 


Input File Usage: 
*ELASTIC, TYPE=ENGINEERING CONSTANTS 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Elasticity > Elastic: Type: Engineering Constants 
Stability 


Material stability requires 
Fy, Fa, Es, Gia, Giz, Go3 > 0, 
Yi2| < (Ei /E2)"”, 
v13| < (E1/E3)"?, 
\v23| < (E2/E3)"”, 


1 — Yy2QV21 — 123132 — 31113 — 22132113 > 0. 


When the left-hand side of the inequality approaches zero, the material exhibits incompressible behavior. Using the relations 
Yj /E,=V;/E;, the second, third, and fourth restrictions in the above set can also be expressed as 


ver] < (E2/E1)"”, 
vst] < (E3/E1)"?, 
|V39| < (sj tay. 


Defining transversely isotropic elasticity 


A special subclass of orthotropy is transverse isotropy, which is characterized by a plane of isotropy at every point in the 
material. Assuming the 1-2 plane to be the plane of isotropy at every point, transverse isotropy requires that #, =F, =Ep, 





V31 =V32="4p, V13=V23=Vpt, and G'13=G23=G:, where p and t stand for “in-plane” and “transverse,” respectively. Thus, while 
tp has the physical interpretation of the Poisson's ratio that characterizes the strain in the plane of isotropy resulting from 
stress normal to it, %¢ characterizes the transverse strain in the direction normal to the plane of isotropy resulting from stress 
in the plane of isotropy. In general, the quantities 44, and Yp_ are not equal and are related by Vtp/ E4=Vpt / Ep. The stress- 
strain laws reduce to 
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eu 1/Ep —p/Ep —Vp/Ex 0 0 0 Cis 
€22 —Vp/Ep 1/Ep —Vp/Ex 0 0 0 On9 
633 | | —vpt/Ep —Upt/Ep 1/E, 0 0 0 O33 
mf | oO 0 0 1G 0 0 on [- 
‘yi3 0 0 0 0 1/G 0 713 
‘23 0 0 0 0 0 1/G; O23 


where Gp=E,/2 (1+ Vp) and the total number of independent constants is only five. 


In Abaqus/Standard spatially varying transverse isotropic elastic behavior can be defined for homogeneous solid continuum 
elements by using a distribution (Distribution definition). The distribution must include default values for the elastic moduli 

and Poisson's ratio. If a distribution is used, no dependencies on temperature and/or field variables for the elastic constants 
can be defined. 


Input File Usage: 


*ELASTIC, TYPE=ENGINEERING CONSTANTS 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Elasticity > Elastic: Type: Engineering Constants 
Stability 


In the transversely isotropic case the stability relations for orthotropic elasticity simplify to 


5, BiG, Gi > 0, 
| <1, 

Vpe | < (Ep/ Ex)”, 
Vep| < (E,/ Ep)”, 


1 — uy? — Wippt — 2WpepVpt > 0. 


Defining orthotropic elasticity in plane stress 


Under plane stress conditions, such as in a shell element, only the values of Ay, Fo, 42, Gi, G43, and G3 are required to 
define an orthotropic material. (In all of the plane stress elements in Abaqus the (1, 2) surface is the surface of plane stress, 
so that the plane stress condition is 733 = 0.) The shear moduli G13 and Gy3 are included because they may be required for 
modeling transverse shear deformation in a shell. The Poisson's ratio 91 is implicitly given as ¥91 = (2/1) V12. In this 
case the stress-strain relations for the in-plane components of the stress and strain are of the form 


E] 1/E, —vp/E, 0 Cu 
€2 => —12/ Ey 1/E» 0 022 
12 0 0 1/Gi2 T12 


In Abaqus/Standard spatially varying plane stress orthotropic elastic behavior can be defined for homogeneous solid 
continuum elements by using a distribution (Distribution definition). The distribution must include default values for the elastic 
moduli and Poisson's ratio. If a distribution is used, no dependencies on temperature and/or field variables for the elastic 
constants can be defined. 


Input File Usage: 


*ELASTIC, TYPE=LAMINA 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Elasticity > Elastic: Type: Lamina 
Stability 


Material stability for plane stress requires 
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Ey, E2, Giz, Giz, G23 > 0, 


\v12| < (Ey /E2)'/?. 


F 


Defining orthotropic elasticity by specifying the terms in the elastic stiffness 
matrix 


Linear elasticity in an orthotropic material can also be defined by giving the nine independent elastic stiffness parameters, as 
functions of temperature and other predefined fields, if necessary. In this case the stress-strain relations are of the form 


(ou ) Din Diz Diis3—O 0 0 E11 Ell 
O22 D2222 D233 0 0 0 E22 E22 
033 D3333 (0 0 0 £33 £33 
012 7 Di22 =O 0 12 7 a 12 
O13 sym Di313. (0 713 713 
023 D323 23 23 


For an orthotropic material the engineering constants define the D matrix as 


Dir = E; (1 — 23132) T, 


( 
Do222 = E> (1 — 4331) T, 
D3333 = E3 (1 — 22) T, 
Dyj22 = Ey (va + 13123) Y = Ee (viz + 3213) TY, 
Dy133 = Fy (v31 + ¥21¥32) TY = E3 (13 + 12023) TY, 
Do233 = Ep (v32 + 1231) T = Es (23 + v2iM3) TY, 
Di212 = Gia, 
Dy313 = Gis, 
D323 = Go3, 
where 
= 1 


1 — 4221 — 493132 — 131113 — 2021132113 


When the material stiffness parameters (the Dijk) are given directly, Abaqus imposes the constraint 033 = 0 for the plane 
stress case to reduce the material's stiffness matrix as required. 


In Abaqus/Standard spatially varying orthotropic elastic behavior can be defined for homogeneous solid continuum elements 
by using a distribution (Distribution definition). The distribution must include default values for the elastic moduli. If a 
distribution is used, no dependencies on temperature and/or field variables for the elastic constants can be defined. 


Input File Usage: 
*ELASTIC, TYPE=ORTHOTROPIC 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Elasticity > Elastic: Type: Orthotropic 
Stability 


The restrictions on the elastic constants due to material stability are 
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Dy111, D2222, D3333, D212, D313, D323 > 0, 


|Dii22| < (Di Deny”, 
|D1133| < (D1111D3333)"””, 
|D2233| < (Do222D3333)"””, 
det (D“) > 0. 


The last relation leads to 


2 2 2 
Dy4111-D2222D 3333 + 2D1122D 11332233 — Do222D — Dy111 D933" — D3333D 499° > 9- 


1133 


These restrictions in terms of the elastic stiffness parameters are equivalent to the restrictions in terms of the “engineering 
constants.” Incompressible behavior results when the left-hand side of the inequality approaches zero. 
ry 
Defining fully anisotropic elasticity 


For fully anisotropic elasticity 21 independent elastic stiffness parameters are needed. The stress-strain relations are as 


follows: 
ou Din Dii22 Dii33 Dire Ditizs D123 E11 E11 
022 D222 D233 D2212 D2213 D2223 E22 E22 
o33 | D3333 D3312 D3313 D3323 33 | (D*} E33 
O12 7 Dj212_ D213 D223 m2 { 12 
013 sym Dy313 D1323 713 713 
023 D323 Y23 ¥23 


When the material stiffness parameters (the Dijk) are given directly, Abaqus imposes the constraint 033 = 0 for the plane 
stress case to reduce the material's stiffness matrix as required. 


In Abaqus/Standard spatially varying anisotropic elastic behavior can be defined for homogeneous solid continuum elements 
by using a distribution (Distribution definition). The distribution must include default values for the elastic moduli. If a 
distribution is used, no dependencies on temperature and/or field variables for the elastic constants can be defined. 


Input File Usage: 
*ELASTIC, TYPE=ANISOTROPIC 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Elasticity > Elastic: Type: Anisotropic 
Stability 


The restrictions imposed upon the elastic constants by stability requirements are too complex to express in terms of simple 
equations. However, the requirement that D? is positive definite requires that all of the eigenvalues of the elasticity matrix 
[DF be positive. 


rt 


Defining orthotropic elasticity for warping elements 


For two-dimensional meshed models of solid cross-section Timoshenko beam elements modeled with warping elements (see 
Meshed beam cross-sections), Abaqus offers a linear elastic material definition that can have two different shear moduli in the 
user-specified material directions. In the user-specified directions the stress-strain relations are as follows: 


oO E E 
T1 = Gi 71 
T2 Go 2 
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A local orientation is used to define the angle @ between the global directions and the user-specified material directions. In 
the cross-section directions the stress-strain relations are as follows: 


- E 0 0 j 
1lp= G1 (cos a)? + Go(sin a)? (G, — G2) cosasina Yi ps, 
T2 sym Gi (sin a)” + G2(cos a)? Y2 


where o represents the beam's axial stress and 7; and 72 represent two shear stresses. 
Input File Usage: 
*ELASTIC, TYPE=TRACTION 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Elasticity > Elastic: Type: Traction 
Stability 
The stability criterion requires that # > 0, G, > 0, and Gy > 0. 


rt 


Defining elasticity in terms of tractions and separations for cohesive elements 


For cohesive elements used to model bonded interfaces (See Defining the constitutive response of cohesive elements using a 
traction-separation description) Abaqus offers an elasticity definition that can be written directly in terms of the nominal 
tractions and the nominal strains. Both uncoupled and coupled behaviors are supported. For uncoupled behavior each traction 
component depends only on its conjugate nominal strain, while for coupled behavior the response is more general (as shown 
below). In the local element directions the stress-strain relations for uncoupled behavior are as follows: 


te == E,, Es 
ty Ex Et 


The quantities t,, t;, and t; represent the nominal tractions in the normal and the two local shear directions, respectively; 
while the quantities €;), €s, and €; represent the corresponding nominal strains. For coupled traction separation behavior the 
stress-strain relations are as follows: 


tn Ean Ens Ent En 
ts = Ens Ess Est Es 
th En Es Ex Et 


Input File Usage: 
Use the following option to define uncoupled elastic behavior for cohesive elements: 
*ELASTIC, TYPE=TRACTION 


Use the following option to define coupled elastic behavior for cohesive elements: 





*ELASTIC, TYPE=COUPLED TRACTION 


Abaqus/CAE Usage: 
Use the following option to define uncoupled elastic behavior for cohesive elements: 


Property module: material editor: Mechanical > Elasticity > Elastic: Type: Traction 

Use the following option to define coupled elastic behavior for cohesive elements: 

Property module: material editor: Mechanical > Elasticity > Elastic: Type: Coupled Traction 
Stability 


The stability criterion for uncoupled behavior requires that En, > 0, Fs, > 0, and Ez > 0. For coupled behavior the 
stability criterion requires that: 


En >90, Ey, >90, Ey > 0; 
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Ens < V Enn Ess} 
Est < V Ess Ett} 
Ent < V Ean Ett} 


Pig, Bs Ent 
det | En, Ess Eg: | > 0. 
Ent Est Ex 


rt 


Defining isotropic shear elasticity for equations of state in Abaqus/Explicit 


Abaqus/Explicit allows you to define isotropic shear elasticity to describe the deviatoric response of materials whose 
volumetric response is governed by an equation of state (Elastic shear behavior). In this case the deviatoric stress-strain 
relationship is given by 


S = 2ue, 


where & is the deviatoric stress and ef! is the deviatoric elastic strain. You must provide the elastic shear modulus, jz, when 
you define the elastic deviatoric behavior. 


Input File Usage: 


*ELASTIC, TYPE=SHEAR 
Abaqus/CAE Usage: 


Property module: material editor: Mechanical > Elasticity > Elastic: Type: Shear 


ry 


Defining the elastic transverse shear modulus 


Abaqus allows you to specify the initial elastic transverse shear moduli for computing the transverse shear stiffness of shells 
and beams. You can define isotropic, orthotropic, and anisotropic transverse shear moduli. 


For beam elements, the material transverse shear stress-relations are as follows: 
e Isotropic: 


\eas=Lo ellos 
0713 0G 713 
e Orthotropic: 
Pere cal 
013 0 Gist Lvs 
e Anisotropic: 
{ \ = ee a ie} 
013 D213 D41313 713 


For shell elements, the material transverse shear stress-strain relations are as follows: 
e Isotropic: 


ee adem 
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e Orthotropic: 


{2} i 0 | AY 
093 0 G3} | 23 


e Anisotropic: 


i \ - ae pal ee } 

023 Dy323  Do323 “Y23 

You can define the transverse shear moduli for materials for which estimates of the shear moduli are unavailable during the 
preprocessing stage of input; for example, when the material behavior is defined by user subroutines UMAT, UHYPEL, UHYPER, 


or VUMAT. You can also define the transverse shear moduli to override the default values computed by a material behavior; for 
example, you can define anisotropic transverse shear moduli for an isotropic elastic material. 


The transverse shear moduli that you define are used only to compute the transverse shear stiffness at the beginning of the 
analysis. Any change to the transverse shear moduli during the analysis is ignored. The transverse shear moduli of a material 
are ignored if nondefault values of the transverse shear stiffness are specified for a section that is associated with the material 
(see Defining the transverse shear stiffness for shell sections, Defining the transverse shear stiffness for general shell sections, 
and Transverse shear stiffness definition for beams). 


Input File Usage: 
Use the following option to specify an isotropic transverse shear modulus: 
*TRANSVERSE SHEAR, TYPE=ISOTROPIC (default) 
Use the following option to specify orthotropic transverse shear moduli: 
*TRANSVERSE SHEAR, TYPE=ORTHOTROPIC 


Use the following option to specify anisotropic transverse shear moduli: 





*TRANSVERSE SHEAR, TYPE=ANISOTROPIC 
Abaqus/CAE Usage: 
Specifying the elastic transverse shear modulus is not supported in Abaqus/CAE. 


rt 


Elements 


Linear elasticity can be used with any stress/displacement element or coupled temperature-displacement element in Abaqus. 
The exceptions are traction elasticity, which can be used only with warping elements and cohesive elements; coupled traction 
elasticity, which can be used only with cohesive elements; shear elasticity, which can be used only with solid (continuum) 
elements except plane stress elements; and, in Abaqus/Explicit, anisotropic elasticity, which is not supported for truss, rebar, 
pipe, and beam elements. 


If the material is (almost) incompressible (Poisson's ratio v > 0.49 for isotropic elasticity), hybrid elements should be used in 
Abaqus/Standard. Compressible anisotropic elasticity should not be used with second-order hybrid continuum elements: 
inaccurate results and/or convergence problems may occur. 
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